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Abstract 

This paper studied the cluster synchronization of directed complex networks with time delays. It is different from undirected 
networks, the coupling configuration matrix of directed networks cannot be assumed as symmetric or irreducible. In order 
to achieve cluster synchronization, this paper uses an adaptive controller on each node and an adaptive feedback strategy 
on the nodes which in-degree is zero. Numerical example is provided to show the effectiveness of main theory. This method 
is also effective when the number of clusters is unknown. Thus, it can be used in the community recognizing of directed 
complex networks. 
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Introduction 

During last decade, the study of complex networks has become 
a hot topic in various fields like physics, mathematics, biology, 
social sciences, computer sciences, and so on [1-3]. Most of 
complex networks have two properties: small-world and scale-free 
[4—5]. Recently, as one of the most important phenomenon of 
dynamical system, synchronization has gained growing attention. 
So far, many different kinds of synchronization in complete 
networks are realized, such as generalized synchronization, phase 
synchronization, cluster synchronization and so on [6-25]. 
Nowadays, cluster synchronization has been widely and thor- 
oughly studied because it can show the community of the complex 
networks [26-30]. 

Cluster synchronization is a middle state of the progress which is 
from none-synchronization to complete synchronization. When 
this middle state is achieved, the nodes in same group (or 
community, or cluster) can achieve complete synchronization, but 
the nodes in different clusters are chaotic. Owing to the significant 
application in biological science and communication engineering, 
the researching of cluster synchronization is focus on the control 
method such as pinning control, adaptive control, impulsive 
control, and so on, but few of them studies cluster synchronization 
of directed complex networks with time delays. 

Liu and others researched generalized synchronization of three 
typical complex dynamical networks including scale-free network, 
small-world network, and a family of interpolating network [7] . 
They found that there is a general progress to global generalized 
synchronization (GS): non-GS —> partial GS — > global GS and the 
GS stats from a small part of hub nodes with larger degrees first. In 
their paper, the partial GS is called cluster synchronization. 
Several interesting adaptive and impulsive synchronization criteria 
are attained for a general complex dynamical network with two 



different clusters by Shi and others in [8]. Lu proposed a novel 
adaptive strategy to make a network achieve cluster synchroniza- 
tion in [9], and Liu and others investigated the cluster 
synchronization with intermittent control in [10]. They also 
pointed out that to realize cluster synchronization, enlarging the 
couplings of nodes in the same cluster is the key point. 

There are also some papers on cluster synchronization of 
directed networks without time delays. Ma and others intensively 
studied the pinning cluster synchronization of directed complex 
networks in [11]. They gave the pinning controllers which are 
applied to inter-act nodes and intra-act nodes with zero in-degree, 
respectively. 

This paper uses an adaptive controller to make a directed 
network with time delays achieved cluster synchronization. The 
rest part of the paper is shown as following. In Section 2, the 
model of directed complex dynamical network and some 
preliminaries are given. The main theorems and corollaries for 
cluster synchronization through adaptive control are given in 
Section 3. At last, a numerical simulation is provided to show the 
effectiveness of the theoretical results. Conclusions are finally 
drawn in Section 5. 

Preliminaries 

Consider a directed complex network with jV identical coupled 
nodes: 

N 

Xi(t)=f i (x i (t),Xi(t-T ] ))+J2 a v x j('- T 2); i=l,2,...,N (1) 

j=i 

Here Xi= [xa,xa,—,Xw] T eR N is the state vector of node z; 
function / :R N ^R N is a nonlinear function which can describe 
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each node's dynamics; Ti and T2 are time-varying delay and 
coupling delay, respectively. Matrix A = (ay) N x N represents the 
topological structure of the network. In a directed complex 
network, ay is defined as follows: if there is a direct link from node 
i to node j(i j), then ay = — = 1 ; otherwise ay = aji = 0. Matrix 
A is satisfied with diffusive condition as follows: 



au=- ^2 a ii> i=l,2,-,N 

Because the complex network is directed, this paper doesn't 
assume A as symmetric or irreducible like other papers. The in- 
degree of node i is defined as: 



Xi(t) =fi(Xi(t),Xj(t - Tl)) + 
N 

aijXjjt — Z2) + Uj(t); i=\,2,...,N 
The controller is designed as following: 



= - X] a 'J Sa >j ( ? ~ T2 ) ~ ke ' 
J=l 



(3) 



(4) 



Here, the constant k>0. The error system of Eq. (2) can be 
obtained as 



-«g= X a y 

7=1 

If the rth node is satisfied with the following equation 



^ay= -a u = 0, 
J=i 

this paper will call the node which is under this condition as 0-in- 
degree node. 

Assume that the network has P clusters (P is unknown), which 
means that all nodes in network will split into P groups when the 
network achieves cluster synchronization. If P=l, then cluster 
synchronization turns to complete synchronization. If node i 
belongs to the kth cluster, this paper denotes that co, = k. When the 
complex network achieves cluster synchronization, for any node 
ij, the following equation is established. 

{lim j j Xi(t) — Xj(t) j j = 0, to, = Wj 
lim||x,-(0-*/(0| #0, coi^coj 

That is to say, when the network achieves cluster synchronization, 
the community of the network can be recognized. If it defines a 
solution vector S = (si(t),S2(f),...,SM(t)) T to represent the desired 
state when the network achieves cluster synchronization at time t, 
here k= 1,2,. ..,M. The error system is defined as follows: 

e i (t) = x i (i)-s mj (t) 

Here the stable dynamic status s w . is satisfied with 
Sa)i(t)=f(s a ,.(t),s a) .(t — T;{)). The complex network can be consid- 
ered to achieve P-cluster synchronization when the following 
condition is satisfied: 



lim 1 1 <?, ( 0 I =0, (Oi = a)j 
lim||e,(0|| #0, coi^coj 



e, = Xj - s ai =f i (x i {t),x i {t - ti )) - 

N ( 5 ) 
fi(So H (t),s m .{t--cx))+ X a i] Xj{t-% 2 ) + Ui 

y=i 

Throughout this paper, the following assumptions are needed to 
prove the main theorem. 

Assumption 1. If there is a nonlinear dynamical function /, 
to any state vectors x,yeR N x 1 , there exists a constant M > 0 to 
make the following equation established: 



\\f(x(t),x(t - 1)) -f(y(t),y(t - 1)) 1 1 < m 
(||x(0-X0ll + ll*O-*)-X'-*)ll) 



8f 



(6) 



Remark 1. Assumption 1 holds as long as — are uniformly 
bounded. Almost all well-known dynamical chaotic and hyper 
chaotic systems have the form of Eq. (3), which meets the 
condition of assumption 1 [31]. 

Assumption 2. There exists a constant fi which can make a 
differentiable time-varying delay x{t) satisfied the following 
equation. 

0<T(0<jU< 1 

It is clearly that assumption 2 is valid for constant z(t). 

1. When the network has no 0-in-degree node. 

Theorem 1. Under assumption 1 and 2, the controlled 
complex network Eq. (3) with adaptive controller Eq. (4) can 
achieve cluster synchronization if k is satisfied with the following 
equation 



(2) 



M+l M+l , , ,- 

Here, M and ji are positive constants, and 
- A A + A T 



(7) 



Adaptive Cluster Synchronization 

In order to make complex network Eq. (1) achieved cluster 
synchronization, an adaptive controller is added on each 

node. The controlled dynamic network can be rewritten as 
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Thus Eq. (2) can be satisfied by the condition as Eq. (7), the proof 
is completed. 

2. When the network has some O-in-degree nodes. 

O-in-degree nodes just send information into network but do not 
receive information from other nodes. That is to say, O-in-degree 
nodes are hard to achieve synchronization. This can be proved as 
following. The linearization of Eq. (5) with controller Eq. (4) can 
be rewritten as 



N N N 

= ^e i T \f i (x i (t),x i (t-T i ))-f i (s 0)i (t),s ai (t-Ti))+^a i jXj(t--z 2 )- ^a ij s aj {t-x 2 )-ke l ] 
1=1 j=\ ,=\ 

M+l T -M(l-ii)^ T , , , , -(1-T 2 )v^ r, 

N M+l M N N N 

< C 9 — eje t + —eJ{t-Ti)e i {t-T X )]+Y j e i T ^ j a ij x j {t-T 2 )-Ya ij s m .(t-T 1 )-ke l ] 

i=i 1 1 i=\ j=\ j=\ 

M+l " T -M(I-ti)^ Tl . , , _(i_i 2 ) * 
M+l M+l f * r * * r , ^ 

y r> i=l i=\ j=\ 

M(ii-fi)^ T -(1-t 2 )v^ T( . . , . 

jv 

e = {Dfi(xi) — ^ a,y)e,- (9) 



Here, Dfi(xi) is the Jacobi matrix of function f(xi) at x,-. If the rth 
node is O-in-degree node, one has 

N 

So O-in-degree node has 

e i = (Df i {x i ))e i . 

Since -D/i(x,)^0, e/T^O, the error system of O-in-degree node is 
hardly to equal 0. It means that the network is hardly to achieve 
synchronization. 

In order to make the network which has O-in-degree nodes 
achieved cluster synchronization, this paper designs a feedback 
adaptive strategy on O-in-degree nodes. 

Theorem 2. Under assumption 1 and 2, if complex network 
Eq. (3) has O-in-degree nodes, it can achieve the desired cluster 
synchronization if controller Eq. (4) and adaptive condition Eq. (7) 
hold, and adaptive feedback strategy is given as Eq. (10), and Eq. 
(11) is established. 



Proof. Define a Lyapunov function: 



2V(t)=J2 e f(t)e,(t) + 



M 



y^e T l {z)e i {z)dz 



+ 



rt j N 



(8) 



ef(z)e,(z)cfe 



Calculating the time derivative of V{i) along the trajectories of Eq. 
(4), one has 



H4 i i 



If the equation E = (ef ,...,e^/) T eR N x 1 is denoted, one gets 



+ W^) -^ T ^AE( t -r 2) + 
l -E T AA T E + 

-(l-t2) £r ^ 

2(1 

According assumption 1, one gets 

E T AE(t-T 2 )< l -E T AA T E+ l -E T (t- x 2 )E{t - x 2 ) 

Thus, under assumption 2 and Eq. (7), the following equation can 
be established. 

M+l M+l T 1 t t 

v ^— + w^r k)EE+ 2 EAAE 

2(1-,) 2(1-,) 

<( M ^ + ^\-k+A m AA)) E r E 
2 2(1-,) 

<0 
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k = m E i 



M+l M+l 



2 2(1- fi) 



-k 0 + A m:lx (A)<0 



(10) 



(11) 



Simulation 

This section will give some examples to verify the effectiveness 
of the proposed theorems in section 3. In the following numerical 
simulations, 3-dimensonal Lorenz system is designed as the 
dynamical of each node. Lorenz function can be described as 
following: 



where m is a positive constant, k§ > 0 is a known constant. 
Proof. Define a Lyapunov function as 



a(x a -x a ) 
/(*;) ={ P x n- x i2 - x n x a 
xnxa-yxa 



(12) 



2V{i)=Y j eJ{i)e i {t)+-{k-k 0 f 

J >-n 1 ^,=i 
(t , N 



+ 



Calculating the time derivative of V{t) along the trajectories of Eq. 
(4), one has 



When the parameters are chosen as a=\0, /? = 28, y = 8/3, 
Lorenz system is chaotic. Under these parameters, the nodes' 
dynamics can be described as 

Xi =fi(xi) = [lOOa - x,i),28x,i —xa -xnxa,x n xa - %/7>xa] T {\7>) 



Example 1 

In this simulation, a directed BA scale-free network is 
constructed. The detail generation algorithm for BA scale-free 
network is introduced in [5]. The parameter are m = mo = 5, 



2(1-/4 jz 



N N N N 

= ^ j ei r [F{xi,t)-F(s m .,t)+^a ij Xj(t-T2)- E ayS^ (f-r 2 )- Are,] +(k-k 0 ) E e, r e,- 

' 7=1 7 = 1 



M+l ^ T 

+ 2(W0§ e ' e ' + - 



2(1-/1) 



i=i 



JV M+1 M JV JV at 

< E [ 7 ef e '+ ^ g f( f - T i) e '( ? - T i)]+ E e ' r E a ;/ x /( f - T 2)- ^fly^^-^)-^,'] 

i=l i=l 7=1 7 = 1 



M+l ^ r 



M(l-i Q^ T(t W( ^ -(l-i 2 )^ r 
2(l-|i)£f-'-* 1 2(1-/1) ._, 



E - Ti)e,0 - • ri) + 2(i £ ^ " T2)e,(t ~ Xl) 



M+l M+l , * T * * r , 

= ( — ^— + 27j—) ~ E e < e i+LL e i <W ~ T2 ) 



iV= 100. Because the network is directed, when the ?th node and 

the jth node are connected from node i to node j, then 

ay = — Oji = 1 . Each node of the network is controlled as Eq. (4), 
According the proof of theorem 1, under assumption 1 and 2, . ■ , r u j ■ t n u- u . u r i, 

& r ' f > anc j m e m-degree ol each node is not 0, which means the following 

one gets V<0 if the following equation is established. equation will be established for each node: 



M+l M+l , . - „ 

— ^— + r - h + A max (A) < 0 

2 2(1 -n) 



100 
7=1 



The proof is complete. 
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Figure 1. The average value of error system in BA scale-free network without 0 in-degree nodes. 

doi:10.1371/journal.pone.0095505.g001 



The attractor of Lorenz system is bounded by |jc,-i|<29, 
|x a |<29, -1<a' /3 <57, |si|<29, \s 2 \<29, -1<j 3 <57, thus 
the network has three clusters. In the following simulation, this 
paper will use the method in section 3 to confirm the number of 
clusters is three. According to theorem 1 , the dynamic network Eq. 
(2) with controller Eq. (4) can achieve cluster synchronization when 



Eq. (7) is established. In order to measure the quality of the process 
of cluster synchronization, this paper uses the following quantities: 



E k (t)= ll*i(0-tt(O|| 2 ; i=h2,-,N; fc = l,2,3. 




15 
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Figure 2. The value of each cluster's stable state in BA scale-free network without 0 in-degree nodes. 

doi:1 0.1 371 /journal.pone.0095505.g002 
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25 
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10 15 20 25 30 35 40 
t 



Figure 3. The average value of error system without adaptive feedback in NW small-world network with 0 in-degree nodes. 

doi:10.1371/journal.pone.0095505.g003 



Here E^(t) represents the average of each error system. When the 
network achieves cluster synchronization, the following equation 
will be satisfied: 



lim Ek(t)-- 

1^0 



= 0; k= 1,2,3. 



The simulation result is shown in the following. It's easy to see 
that each error system is 0 at last which means that each cluster 
achieves synchronization from Fig. 1. Fig. 2 shows the value of 
Sk(t). Because s\ =£sj ¥=si when E\ =Ej = £3 =0, it is shown that 
the nodes in different cluster cannot achieve synchronization 




i 1 1 1 r 



_i 1 1 i_ 



10 15 



20 
t 




25 30 35 40 



Figure 4. The average value of error system with adaptive feedback in NW small-world network with 0 in-degree nodes. 

doi:10.1371/journal.pone.0095505.g004 
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Figure 5. The value of each cluster's stable state with adaptive feedback in NW small 

doi:10.1371/journal.pone.0095505.g005 



-world network with 0 in-degree nodes. 



clearly. The network has three clusters according to the simulation 
result. 

Example 2 

In this example, this paper uses a directed WS small-world 
network. The number of nodes as the BA scale-free in example 1 . 
In detail, this paper will use the parameter to construct a WS 
small-world network as [4], the rewiring probability is /> = 0.2, the 
number of nodes is 100, and k = 5. In this example, there are some 
nodes will be chosen randomly as 0-in-degree nodes. 

At first, controller Eq. (4) is added on each node. The simulation 
result is shown as Fig. 3. It is easy to see that each cluster can't 
achieve synchronization at all because each error system cannot 
achieve 0. Then the adaptive feedback strategy Eq. (10) is added 
on each 0-in-degree node, here. The simulation result shows that 
the network can achieve cluster synchronization as Fig. 4 and 
Fig. 5. It is easy to see that the network can achieve cluster 
synchronization when the adaptive feedback strategy as Eq. (10) is 
used. It is easy to see that the number of clusters is three in Fig. 5. 

Conclusion 

In this paper, cluster synchronization of directed complex 
dynamic network with time delays was investigated. An adaptive 
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